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D ' Received 

Aims. Gamma-ray line emission from the radioactive decay of '^''Al reflects nucleosynthesis in massive stars and super- 
[j^ [ novae. We use INTEGRAL ^®A1 measurements to characterize the distribution and characteristics of "^"Al source regions 

p3 , throughout the Galaxy. 

HH ' Methods. The spectrometer SPI aboard INTEGRAL has accumulated over five years of data on ^®A1 gamma-ray emis- 

(-H ' sion from the Galactic plane. We analyzed these data using suitable instrumental-background models and adopted sky 

distribution models to produce high-resolution ^®A1 spectra of Galactic emission, spatially resolved along the Galaxy 
plane. 

Q ■ Results. We detect the ^''Al line from the inner Galaxy at ~ 28a significance. The line appears narrow, and we con- 

V-^ ' strain broadening in the source regions to < 1.3 keV (2a"). Different sky distribution models do not significantly affect 

C/j I those large-scale results. The ^^Al intensity for the inner Galaxy is derived as (2.9 ± 0.2) x 10~* ph cm""^ rad~^, 

, consistent with earlier results from COMPTEL and SPI data. This can be translated to an ^^Al mass of 2.7 ± 0.7 

Mq in the Galaxy as a whole. The ^®A1 intensity is also confirmed to be somewhat brighter in the 4th than in the 
1st quadrant (ratio ~ 1.3 ± 0.2). ^®A1 spectra separately derived for regions along the Galactic plane show clear line 
centroid shifts, attributed largely to the Galaxy's large-scale rotation. The ^®A1 line toward the direction of the Aquila 
region (20° < / < 40°) appears somewhat broadened. Latitudinal variations of "^^Al emission towards the inner Galaxy 
are studied, finding a latitudinal scale height of I30I7Q'' P*^ (!"") ^'^^ ^^Al in the inner Galaxy and a hint (3(t) of peculiar 
, ^®A1 emission towards the region I < 0° , 6 > 5°. 

, Key words. Galaxy: abundances - ISM: abundances - nucleosynthesis - gamma-rays: observations 
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^ 1. Introduction The detailed study of ^^Alline emission from the Galaxy 



is one of the main science goals of the INTEGRAL mission. 
The unstable isotope '^^Al has a mean hfetime of 1.04 Myr. SPI aboard INTEGRAL is a high-resolution spectrometer 
26A1 first decays^ into an excited state of '^^Mg, which de- with energy resolution of 3 keV (FWHM) at 1809 keV, 
excites into the ^'^Mg ground state by emitting gamma-ray which therefore adds high-resolution spectroscopic infor- 



O 

> 



H , photons with the characteristic energy of 1809 keV . The mation to ^^A\ astronomy The Compton Observatory in 

■ - - • 1809 keV gamma-ray line from radioactive ^^Al serves as a its early design phase included a fifth instrument 'GRSE', 

tracer of the recent nucleosynthesis activity in the Galaxy a high-resolution spectrometer later abandoned for cost and 

The 1809 keV 7-ray line emission from the Galaxy was complexity reasons. COMPTEL as a scintillation-detector 

first detected with the Ge spectrometer on the HEAO-C based instrument had modest spectral resolution of about 

space experiment (Mahoney et al. 1982). The Compton 10% (FWHM). The detailed measurement of ^""Al hne po- 

Observatory sky survey 1991-2000 then with COMPTEL sition and shape is expected to reveal more information be- 

imaging of the ^®A1 line across the sky showed that ^^Al yond the COMPTEL imaging survey about the ^^Al sources 

emission extends along the Galactic plane, thus clearly and their location through the Doppler effect, induced from 

establishing ^^Al nucleosynthesis as a widely-distributed Galactic rotation and dynamics of the ejected ^^Al as it 

Galactic phenomenon (Diehl et al. 1995, Pliischke et al. propagates in the interstellar medium around its stellar 

2001). The structure of this emission, alignments of emis- sources, 
sion maxima with spiral- arm tangent, and comparisons 

with tracers of candidate ^SAl sources, all point to a pre- Earher INTEGRAL analysis had used 1.5 years of SPI 

dominant origin of ^^A\ in massive stars (Prantzos and data to first explore the large-scale spectral characteristics 

Diehl 1996, Chen et al. 1995, Diehl et al. 1995, Knodlseder of Al emission in the inner Galaxy (Diehl et al. 2006a, 

2999^ 2006b). A detection of the Al line from the inner Galaxy 

with a significance of ~ 16a confirmed the narrowness of 

the 26A1 hne (FWHM < 2.8 keV, 2a), which had already 

Send offprint requests to: W. Wang, wwangQmpe.mpg.de been seen by RHESSI (Smith et al. 2003) and HEAO-C 
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Fig. 1. Exposure of the sky in Galactic coordinates (the 
number at the color bar in units of ksec) for the data 
selected from 5-year SPI observations for our ^^Al study 
(INTEGRAL orbits 43 - 650). The database covers the 
whole sky, with 29736 individual pointings, equivalent 
to a total (deadtime-corrected) exposure time of 61 Ms. 
Observations emphasize the Galactic plane, specially in the 
inner Galaxy region, but also specific regions of interest 
such as Gas A, Cygnus, Carina- Vela, Orion, Virgo, and the 
Crab pulsar and nebula. 
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Fig. 2. The unfihered total GEDSAT rates and background 
rates for detector around 1809 keV in model fittings versus 
time (the INTEGRAL Juhan Date, starts at 1 Jan 2000). 
The total GEDSAT rates per detector generally vary in a 
range from 150 - 300 counts/s, but some erratic cases like 
zero values and solar flare spikes are excluded from science 
analysis with our primary selections (also see the text). 



(Mahoney et al. 1984) earlier. We use INTEGRAL/SPI 
data accumulated from five years to extend this study to- 
wards spatially-resolved details of ^^Al line spectroscopy 
across the inner regions of the Galaxy. In this paper, we 
concentrate on the spectral and intensity variations of ^^Al 
emission along the Galactic plane. The main goal is to ex- 
plore if there are variations from the large-scale properties 
of ^^Al emission, such as bulk motion or enhanced turbu- 
lence in specific regions. Our study will be limited by the 
brightness of the signal per region, and hence proceed from 
broad to more constrained regions along the plane, as the 
current exposure and ^^Al brightness allows; with more ex- 
posure, in particular towards outer regions of the Galaxy at 
higher longitudes, the extended INTEGRAL mission is ex- 
pected to eventually allow such studies up to the spatial res- 
olution of about 2.8° (FWHM) of the SPI instrument. We 
also derive spectra for separate regions of Galactic latitudes, 
to probe the symmetry and the scale height of Galactic ^^Al 
emission. 

In this paper, we first describe the SPI observations 
and methods of data analysis. Then we proceed with an 
update of the large-scale characteristics of Galactic ^^Al 
emission, before refining the spatial resolution of our study. 
We discuss the implications of our findings with respect to 
the large-scale distribution of ^^Al sources, the interstellar 
medium in their vicinity, and possible regional deviations. 
Studies of specific regions and their ^^Al emission are un- 
derway and will be reported in separate papers (Cygnus - 
Martin et al., in preparation; Orion - Lang et al., in prepa- 
ration; Sco-Cen - Diehl et al., in preparation). 

2. Observations and Data Analysis 

2.1. SPI observations 

The INTEGRAL spacecraft was launched on October 17, 
2002, into a high-inclination, high-eccentricity orbit in- 
tended to avoid the increased background from the Earth's 



trapped radiation belts. INTEGRAL'S orbital period is ^ 
3 days. The spectrometer SPI consists of 19 Ge detectors 
actively shielded by a BGO anti-coincidence shield. It has 
a tungsten coded mask in its aperture which allows imag- 
ing at ~ 2.8° resolution within a 16° x 16° full coded field 
of view (imaging on INTEGRAL is mainly performed at 
lower energies by the IBIS telescope, with which SPI is 
co-aligned). The Ge detectors are sensitive to gamma-rays 
between 15 keV and 8 MeV, with a total effective area 
^ 70 cm^ at 1 MeV, and achieve an energy resolution of 
- 2.5 keV at 1 MeV (Roques et al. 2003, Attie et al. 2003). 
However, cosmic-ray (CR) impacts degrade this resolution 
over time, and the instrument SPI is periodically switched 
off for 14 days twice a year while annealing (by heating from 
cryogenic temperatures ~ 80 K to 100°C) is applied to the 
detectors to restore the energy resolution back to its pre- 
launch value by thermal curing heating of the CR-induced 
defects (Roques et al. 2003, Leleux et al. 2003). 

In space operations, INTEGRAL with its IBIS and SPI 
telescopes is pointed with a fixed attitude for intervals of 
typically 2000 s (referred to as pointings) , which are suc- 
cessively arranged as a standard pattern of neighbouring 
pointings ~ 2° apart (dithering), and covering target re- 
gion of interest for improved imaging (Jensen et al. 2003; 
Courvoisier et al. 2003). 

We obtain a database from the 5-year SPI observations, 
i.e. from the INTEGRAL orbits 43 - 650, which encom- 
passes 29736 pointings of the spacecraft and its instruments 
across the sky, equivalent to a total deadtime-corrected ex- 
posure time of 61 Ms (see the exposure map in Figure 1). 

The instrument operation has been interrupted from 
a few short anomalies and the occurrences of solar flares, 
and most significantly for the typically two-week annealing 
episodes described above, and by the regular perigee pas- 
sages with switch-offs due to the high background intensity 
from the radiation belts. The sensitivity of our observa- 
tions were further reduced by the failure of two of the 19 
detectors (December 2003 and July 2004). 
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Fig. 3. The COMPTEL maximum entropy ^SAl map (left, Pliischke et al. 2001, hereafter COMPTEL MaxEnt) and a 
homogenous disk model (right, hereafter Homo Disk, constant brightness along longitudes, exponential-like in latitudes 
with scale height 200 pc) for the inner Galaxy as the input sky models in the model fitting. 
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Fig. 4. Residuals versus time (pointings, per 30 min, left) and energies (right) after model fitting of the inner Galaxy region 
for the SE database. Residuals around zero confirm that our background models are adequate, with /d.o.f. = 0.958 in 
model fittings. 



2.2. Data Analysis 

Our data analysis steps include: (1) selection and assem- 
bly of data which are free of contaminations by anomalies 
and, e.g., solar-flare events; (2) modelling the instrumental 
background; (3) deriving spectra by fitting the measured 
data in narrow energy bins with our background model and 
a model of the spatial distribution of celestial gamma-ray 
emission, folded through the instrumental response into the 
data space of the measurement. These steps have been de- 
scribed in several papers, e.g. Strong et al. (2005), Diehl et 
al. (2003), and as applied to gamma-ray lines specifically in 
Wang et al. (2007) and Diehl et al. (2006b). Here we give 
a concise summary and add comments as relevant for the 
studies of this paper. 

We first apply a selection filter to reject corrupted or 
invalid data. This is done for data sections called 'science 
window', typically a time interval of ~ 30 minutes corre- 
sponding to one pointing. We apply selection windows to 
'science housekeeping' parameters such as the count rates in 
several background tracers, proper instrument status codes, 
and orbit phase. For tracing background, we employ the SPI 
plastic scintillator anticoincidence counter (PSAC), and the 
rate of saturating events in SPI's Ge detectors (i.e., events 
depositing > 8 MeV in a single Ge detector; hereafter re- 
ferred to as GEDSAT rates, see Figure 2). This leads to 
exclusions of solar-flare periods and other periods of clearly 
increased / abnormal backgrounds. Additionally, regular 
background increases during and after passages through the 
Earth's radiation belts are eliminated through a 0.05-0.916 
window on orbital phase. 



From the selected events, spectra are accumulated per 
each detector and pointing in energy bins of 0.5 keV in the 
spectral range of 1785 - 1826 keV, and together with dead 
time and pointing information assembled into the analysis 
database. From identically-selected data, we also establish 
a database for the background modelling from the adja- 
cent continuum, combining broader spectral ranges on both 
sides of the ^^Al line into single bins. 

We select events which trigger one and only one of our 
19 Ge detectors ('single events', SE). We exclude the ~50% 
fraction of 'multiple' events, because for events interacting 
in more than one detector different combinations of single- 
detector energies with their correspondingly-different reso- 
lutions would be superimposed for identical values of total 
energy, thus leading to an ill-determined spectral energy 
response. This reduces the total signal obtained, but guar- 
antees that the spectral response of the instrument is well- 
defined for our dataset. 

Instrumental background is dominated by a rather 
smooth spectral continuum, with instrumental lines super- 
imposed - for our case, a broad instrumental-line feature 
centered at 1810 keV comprises 30% of the count rate (see 
Fig. 2 in Diehl et al. 2006b). For our model of the continuum 
background, we make use of the simultaneously-measured 
events in energy bands adjacent to the ^®A1 line region, 
at 1785 - 1802 keV plus 1815 - 1826 keV, and use each 
detector's count rates in this adjacent range to normalize 
the GEDSAT time series with their superior statistical pre- 
cision at regular intervals of several days. In Figure 2 we 
show both the original GEDSAT count rate history and the 
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background model for one of our detectors, as it varies over 
the dataset. It is seen that anomalous spikes of our back- 
ground are excluded from our data, and that there still 
is significant temporal structure in our fitted background 
model. This 'adjacent-energies' background variation model 
per detector thus is based on high-statistics GEDSAT rates 
over short time scales of ~ 100 pointings, as these have 
been found to trace the prompt CR activation rather well. 
Adjusting these 'templates' to the actual counts per detec- 
tor in the 0.5 keV bins used for spectral analysis of the sky 
signal then ensures that over longer time scales (typically 3 
days, where counts in those narrow bins are statistically suf- 
ficient) any second-order deviations from GEDSAT tracing 
actual backgrounds are accounted for. Through this pro- 
cedure we ensure that pointing-to-pointing variations are 
given by our background tracer, independent of the coded- 
mask orientation on the sky. The broad instrumental line 
feature centered at 1810 keV is actually found to also be 
rather well modelled by these 0.5 keV bin adjustments, indi- 
cating that the line and continuum backgrounds have sim- 
ilar temporal behaviour. 

The cosmic-signal spectra arc obtained from our mea- 
sured database of detector spectra when we combine above 
background models with a spatial model for sky emission 
(Figure 3), allowing for adjustments of intensity parameters 
for background and sky intensities per energy bin: 

m,n j=l t i=l 

where e,d,p arc indices for data space dimensions: en- 
ergy, detector, pointing; m,n indices for the sky dimensions 
(galactic longitude, latitude); A is the instrument response 
matrix, / is the intensity per pixel on the sky, fci is the num- 
ber of independent sky intensity distribution maps; is the 
number of background components, S is the count residue 
after the fitting. The coefficients for the sky map inten- 
sity are constant in time, while 0b,t is allowed time depen- 
dent, to cater for different background normalizations for 
each camera configuration of 19/18/17 functional detector 
elements. The sky brightness amplitudes Ps comprise the 
resultant spectra of the signal from the sky. We generally 
use a maximum-likelihood fitting method (implemented in 
a SPI standard tool called spim,odfit, properly accoimting 
for Poisson statistics in our spectra; for more details, see 
Strong et al. 2005). 

We thus obtain per energy bin the fitted intensity pa- 
rameter values with their uncertainties, their covariance 
matrices, and the fit residuals. Analysis of residuals, un- 
certainties, and covariance matrices are made to determine 
the validity of the fit of such a combined model to our 
measured set of spectra. Residuals after the model fitting 
are shown in Figure 4 (residuals with the time and ener- 
gies). Reduced values around 1.0 confirm that both our 
background model and the sky distribution model(s) are 
adequate. 

We use the sky intensity distribution of ^^Al as derived 
from 9 years of COMPTEL observations (maximum en- 
tropy map, Pliischke et al. 2001) as the standard reference 
sky model. We compare this to different ^^Al emission tracer 
models or maps, to assess the impact of different sky dis- 
tribution models on the ^^Al line shape and intensity (see 
§3). For simultaneous fits of sets of partial-sky distribution 



models (in order to derive spatially-resolved spectra), we 
use analytical forms of sky distribution models with pre- 
scribed symmetry to avoid biases to such analysis (see be- 
low). 

After we obtain fitted-amplitude coefficient spectra for 
the sky components, we characterize these spectra and in 
particular the ^^Al line details therein through two differ- 
ent approaches. When the sky signal is weak, we fit the 
spectra with Gaussians plus a linear residual background, 
and use the Gaussian intensities, centroid energies, and 
FWHM widths for relative comparisons, such as trends 
along the plane of the Galaxy. When the signal is suffi- 
ciently strong so that we are sensitive to line width details, 
we describe the line component not by a single Gaussian 
any more, but rather by the convolution of a Gaussian 
with the asymmetric instrumental line response as it de- 
velops from degradation of our detectors' resolution and 
their periodic restorations through annealings. This allows 
us to infer immediate information about the celestial ^^Al 
dynamics, which we identify with this Gaussian, and in par- 
ticular its width, which arises from Doppler shifting of the 
line energies with motion of decaying ^''Al nuclei relative to 
the observer. In both cases, we fit these spectral models to 
the amplitude values and their uncertainties as determined 
from the model fit to the large set of spectra, performing 
a maximum-likelihood fit with the Levenberg-Marquardt 
algorithm. In the latter case, we derive the entire proba- 
bility distribution of the fitted parameter values through 
a Monte Carlo Markov Chain (MCMC) method, which al- 
lows us in particular to determine probability constraints 
for the celestial ^^Al line broadening from kinematics. This 
parameter challenges the spectral resolution of our instru- 
ment (about 3 keV at the ^^Al line energy), hence so far 
is bounded from above only. Such asymmetric parameter 
probability distributions may be far from Gaussian, hence 
require such more sophisticated treatment (Kretschmer et 
al. 2009, in preparation). 



3. °AI in the Inner Galaxy - An Update 

The spectral characteristics of ^^Al emission in the inner 
Galaxy serve to study the current nucleosynthesis activity 
and the properties of the interstellar medium near the ^^Al 
sources on a large-scale averaged scale. We define the "inner 
Galaxy" as the region -30° < / < 30°, -10° < 5 < 10°), 
and may use this as a representative region for this purpose, 
since it coincides with the bright ridge of observed 1809 keV 
emission as observed along the plane of the Galaxy. We 
specifically exclude regions at higher longitudes such as 
Cygnus: This region in particular has been recognized as 
special, in its supernova to Wolf-Rayet ratio, for example 
(Knodlseder et al. 2004); other emission at large longitudes 
also may be deviant from the Galactic average and over- 
emphasized because nearby. We therefore believe that using 
the —30° < I < 30° region will give us a more representa- 
tive picture of ^^Al source environments in the Galaxy. 

Fitting our set of observations to the sky intensity dis- 
tribution of the ^^Al maximum entropy image (MaxEnt) 
from COMPTEL (Figure 3, left), we obtain the updated 
inner- Galaxy ^^Al emission spectrum shown in Figure 5. 
MaxEnt is one possible sky intensity distribution compati- 
ble with the COMPTEL data. The ^^Al Une is detected at 
~ 28cr significance. 
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Fig. 5. Spectrum derived from sky model fitting using the COMPTEL ^^Al Maximum Entropy image. The left figure 
shows the ^^Al line fitted with a Gaussian, the width of ~ 3.16 keV (FWHM) being consistent with instrumental line 
widths around 1.8 MeV. The right figure shows the line fitted with a composite line-shape model using the time-averaged 
instrumental response as it results from cosmic-ray degradation and annealings during the time of our measurement, 
convolved with a Gaussian representing the cosmic (intrinsic) ^^Al line width. The latter is found to be (< 1.3 keV, 2a). 
Both fits find that the line is intrinsically narrow. Systematic variations of derived line fluxes using two spatial models 
are less than the statistical uncertainty in the measurement (fluxes are quoted in units of 10~^ ph cm^^ s^^ rad~^) 



The ^^Al gamma-ray flux from the inner Galaxy 
turns out as (2.93 ± 0.15) x 10""* ph cm^^ s~^ rad~^. 
This is consistent with our earlier values (3.3 ± 0.4) x 
lO"'^ ph cm-2 rad^i (Diehl et al. 2006a, 2006b), and 
also with the COMPTEL imaging- analysis value of (2.8 ± 
0.4) X 10"'' ph cm-2 s-i rad^^ (Pliischke et al. 2001). 
The above value is derived using an asymmetric line shape 
as best matching our expectations from SPI's spectral re- 
sponse and eventual additional celestial line broadening 
the total ^^Al gamma-ray flux of the Gaussian fit as de- 
termined for the inner Galaxy region is (2.73 ± 0.17) x 
10^^ ph cm^^ s^^ rad~^. 

In our model-fitting approach to derive ^^Al spectra 
from SPI spectra per pointing (Eq. 1), the adopted sky 
distribution models may affect the ^^Al flux and line shape 
results. We estimate the variations and potential systematic 
uncertainties introduced by the spatial distribution model 
of the ^^Al emission through variation of this sky distribu- 
tion model within a plausible range. 

Our "standard" is the measurement of ^^Al gamma- 
ray emission directly, from COMPTEL. Uncertainties from 
COMPTEL imaging analysis have however been shown 
to allow for a range of images which are consistent with 
the COMPTEL measurements (Knodlseder et al. 1999). 
We consider the MaxEnt image the compromise between 
making use of COMPTEL's imaging resolution of about 
3.8° (FWHM) and the attempts to suppress artifacts 
from statistical noise. A more conservative COMPTEL 
image is obtained from the Multi-resolution Expectation 
Maximization method (MREM) , which carefully eliminates 
noise contributions at each iteration and builds up the im- 
age starting from large spatial scales, terminating once the 
image obtained does not need further refinement in this 
statistical sense (Knodlseder et al. 1999; Knodlseder 1999). 
The smooth MREM image from COMPTEL also shows the 
inner Galactic ridge as well as Cygnus being bright in ^^Al 
emission, yet does not show the few-degree scale features 



of the MaxEnt map; we consider these two maps as ade- 
quate tests for systematics from the range of direct ^^Al 
measurements. 

Considering the limitations of gamma-ray telescopes, 
it has been plausible to alternatively use maps obtained 
in astronomically-more developed wavelength bands, once 
it it clear that those trace ^^Al sources in the Galaxy. 
Detailed studies have shown (Knodlseder et al. 1999; Diehl 
et al. 1996) that among the best tracers of ^^Al sources 
are (1) the infrared emission of warm dust grains, mapped 
with the COBE/DIRBE and arising from radiative heat- 
ing of dust around clusters of massive stars (Bennett et al. 
1996), and (2) radiation of free electrons (free- free emis- 
sion, Bremsstrahlung) observed at radio frequencies with 
the WMAP satellite, arising from the ionizing massive-star 
radiation around massive star clusters. From astrophysical 
arguments, also maps of interstellar gas in different forms 
should trace the locations and space density of ^^Al sources. 
(3) Molecular gas is observed through CO line emission at 
radio frequencies, and has been mapped in rather fine res- 
olution (Dame et al., 1987 and 2001), (4) atomic hydrogen 
(HI) sky surveys have been accumulated (e.g. Dickey and 
Lockman 1990), and (5) cosmic-ray interactions with in- 
terstellar gas produces penetrating continuum gamma-ray 
emission which has been mapped in the EGRET sky survey 
(Hunter et al. 1997). 

Finally, analytical models for the distribution of ^^Al 
sources have been constructed, based on above knowledge 
of Galactic structure in its different components (Robin et 
al. 2003), properly weighted from astrophysical arguments. 
(6) Double-exponential functions (in galactocentric radius, 
and scale height above the Galactic plane) have been con- 
structed, as well as more sophisticated models including (7) 
spiral structure components and building on the distribu- 
tion of free electrons in the Galaxy as derived from pul- 
sar dispersion measurements (Taylor & Cordes 1993, and 
Cordes & Lazio 2002). The scale height of ^^Al sources has 
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been fomid to lie between the molecular disk (about 50 pc) 
and the thick disk (about 0.3-1 kpc), with plausible values 
around 200 pc. 

We compare ^®A1 line spectra determined from these 
different models and tracers of ^®A1 sources in the Galaxy. 
Fifteen different sky distribution maps have been analyzed, 
and variations on ^^Al brightness, ^^Al line centroid and 
width parameters are shown in Figure 6. Systematic vari- 
ations are within statistical uncertainties, when we vary 
the spatial models for ^^Al emission. We use their scatter 
to estimate a "systematic" uncertainty, which turns out as 
(2.9 ± 0.2) X 10-^ph cm-2 s"! rad"!, 1809.0 ± 0.09 keV, 
0.5 ± 0.45 keV for the flux, centroid, and width parameters 
(with l(j error bars), separately. 

As before (Diehl et al., 2006a), we convert our mea- 
sured ^^Al intensity into an estimate of the total cur- 
rent ^^Al mass in the Galaxy, using an assumed geometri- 
cal source-distribution model to extrapolate across the en- 
tire Galaxy from our inner- Galaxy normalization of such 
a model, as discussed above. This is required because the 
flux to mass conversion relies on ^^Al source distances, not 
measured directly in projected sky brightness distribution 
maps. Several three-dimensional distribution models are 
applied and compared here, e.g., an exponential disk model, 
a "young-disk" model, a geometrical representation based 
on dust emission, and a multi-component model including 
spiral-arm structures and based on abundances of free elec- 
trons in interstellar space. We vary scale height parameters 
of the appropriate components to also study the latitude 
extent of the ^^Al emission. In all models, we have taken 
the distance of the Sun to the Galactic center as Rq = 
8.5 kpc (^®A1 mass sensitively depends on Rq, smaller Rq 
will globally reduce the size of the Galaxy, and result in the 
smaller amount of ^^Al ). 

The average measured ^^Al flux of (2.9 ± 0.2) x 
10~* ph cm^^ s^^ rad^^ for the inner Galaxy thus trans- 
lates into a Galactic ^^Al mass of (2.7 ± 0.7) Mq using 
a plausible scale height of 180 pc. If we ignored the ejec- 
tion of ^^Al into surrounding cavities and corresponding 
champagne flows (see our own scale height determination 
below), and used the lower scale heights of O stars or of 
the molecular disk, we would obtain lower total amounts 
around or even below 2 Mq . This emphasizes the need 
for spatially-resolved ^^Al studies (see Sect 5 and Fig. 12 
below). The quoted uncertainty here includes both the sta- 
tistical uncertainty as propagated from the fitting method, 
and the systematical uncertainty derived from variations 
among plausible representations of the sky distribution of 
the emission plus alternative 3-dimensional ^^Al source dis- 
tribution models (see also Figure A4.1 in Supplemental ma- 
terials of Diehl et al. 2006a). 

The line width of ^^Al conveys global information 
about the spread in projected velocities of ^^Al nuclei when 
they decay and emit characteristic 1809 keV gamma-rays. 
Astrophysical origins of ^^Al line width may come from two 
effects: random motions in the interstellar medium (Chen et 
al. 1997) and Galactic differential rotation (Kretschmer et 
al. 2003). A line broadening of ^ 1 keV from our line-shape 
constraints corresponds to thermal Doppler velocities of ^ 
120 km s-^ 

The Gaussian used to represent the ^^Al line (Fig. 5 
left) shows a width of 3.16 ± 0.15 keV (FWHM), consis- 
tent with the instrumental width determined as ^ 3.1 keV 
near 1809 keV from nearby instrumental lines at 1764 and 



1779 keV. This indicates that the cumulative ^^Al line emis- 
sion in the inner Galaxy is intrinsically rather narrow in 
energy around the laboratory value, with little kinematic 
Doppler broadening, well below SPI's energy resolution. 

For such a strong signal, we can be more ambitious, 
however, and attempt to quantitatively constrain the in- 
trinsic (celestial) line broadening in a statistical probability 
analysis as described above (see also Diehl et al. 2006b and 
Kretschmer ct al. 2004). We know the instrumental line re- 
sponse of our instrument at each epoch of the mission from 
detailed studies of a large number of instrumental lines. 
The response gradually degrades and develops an increas- 
ing tail component on the low-energy side of the photopeak 
energy, due to incomplete charge collection as the detector's 
Ge lattice experiences damaging from cosmic-ray impacts. 
Annealings restore better charge collection and a symmetric 
instriuTiental-line response whose width is determined and 
dominated by statistics of the charge collection. We assem- 
ble the effective instrumental-line response for our study 
through appropriate weightings of these time-dependent re- 
sponses. Then we fit the convolution of a parameterized 
Gaussian with this response shape function, and thus derive 
values for the celestial ^®A1 line centroid energy and broad- 
ening. MCMC sampling of the probability distributions of 
these parameters then allows us to translate the parameter 
value uncertainties into probability constraints by integrat- 
ing over these distributions. In particular, the probability 
distribution of the intrinsic width may not peak at zero, 
indicating a small but by itself non-detectable broadening 
of the line. In this case, our approach yields a reliable esti- 
mate of the upper bound on line broadening, which is the 
astrophysically-relevant quantity to constrain kinematics of 
decaying ^^Al nuclei in the interstellar medium. The fitted 
parameters are the line centroid, the intrinsic width of ce- 
lestial ^^Al , the intensity of the line, and two parameters 
for the underlying continuum (see Figure 5 right). The in- 
trinsic line width is now constrained to < 1.3 keV (2 cr, 
which is consistent and yet significantly smaller than the 
earlier constraints based on fewer SPI data (Diehl et al. 
2003, 2006b). 

In the inner Galaxy, Galactic differential rotation alone 
can lead to significant Doppler shifts towards specific longi- 
tudes where the projected- velocity differences with respect 
to the solar orbit reach maxima; line broadening results 
if we integrate over a larger longitude range with different 
bulk velocity differences. Kretschmer et al. (2003) have sim- 
ulated the ^^Al line shape diagnostics in the inner Galaxy 
due to Galactic rotation and ^®A1 ejection from sources, 
and find that line broadening of up to 1 keV is expected if 
the signal is integrated over the inner region of the Galaxy. 
Our present large-scale line-shape constraints are consistent 
with these expectations. 

If we int(Tprc;t line broadening of the ^''Al line from the 
inner Galaxy in terms of interstellar-medium characteris- 
tics, the intrinsic-width constraint of < 1.3 keV corresponds 
to 160 km s~^ as a corresponding 2a limit on ISM veloci- 
ties. This is well within the plausible and acceptable range 
for the environment of normal interstellar-medium turbu- 
lence (Chen et al. 1997), leading us to conclude that, within 
uncertainties, the average velocities of decaying ^^Al in the 
Galaxy are not abnormally-high (compare discussion after 
Naya et al. 1996 in Chen et al. 1997). 

In summary, the measured line width of ^''Al from the in- 
ner Galaxy is consistent with Galactic rotation and modest 
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Fig. 6. Comparison of the three hne parameters (centroid energy, flux, width, with 1 a error bars) derived in model 
fittings for different distribution models: 1. COMPTEL maximum entropy ^^Al emission map (Pliischke et al. 200f); 2. 
COMPTEL MREM ^^Al emission map (Pliischke et al. 200f); 3. an exponential disk model (scale radius 4 kpc, scale 
height 180 pc); 4. a homogenous disk model (Fig. 3); 5. HI (Dickey & Lockman 1990); 6. CO (Dame et al. 1987); 7. radio 
408 MHz (Haslam et al. 1995); 8. DIRBE/COBE 24Qfim (Bennett et al. 1996); 9. IRAS 12/im (Wheelock et al. 1991); 
10. EGRET (> 100 MeV, Hunter et al. 1997); 11. a young disk model (Robin et al. 2003); 12. the free electron density 
distribution model (TC93, scale height 150 pc, Taylor & Cordes 1993) 13. TC93, scale height 300 pc; 14. the free electron 
density distribution model (NE2001, scale height 140 pc, Cordes & Lazio 2002); 15. NE2001, scale height 330 pc. 




1800 1805 1810 1815 1820 1800 1805 1810 1815 1820 

Energy [keV] Energy [keV] 

Fig. 7. ^^Al spectra for two Galactic quadrants (left 0° < Z < 60°, and right —60° < I < 0°). Line centroids relative to 
the centroid energy of ^^Al line in the laboratory (1808.65 keV) show a significant blueshift in the 4th quadrant. Both the 
spectra have width values near the value of the instrumental line width, implying that ^^Al emissions from two quadrants 
are the narrow lines. In addition, the ^^Al flux of the 4th quadrant is higher than that of the 1st quadrant, and the flux 
ratio is ^ 1.3. 



interstellar-medium turbulence around the sources of ^^Al 
. This confirms earlier results on the ^®A1 line width from 
HEAO-C (Mahoney et al. 1984), RHESSI (Smith 2003), 
and SPI on INTEGRAL (Diehl et al. 2006b). The GRIS 
balloon experiment had reported a very broad line with a 
width ~ 5.4 keV (Naya 1996), which is inconsistent with 
these measurements and clearly ruled out. This is reassur- 
ing: within uncertainties, the average velocities of decaying 
^^Al in the Galaxy are not abnormally-high (see discussion 
in Chen et al. 1997). 

The relative position of the ^®AI gamma-ray Hne 

centroid with respect to the laboratory-determined -^"^Al 
decay energy value of 1808.65 ±0.07 keV (Firestone & 
Ekstrom 2004) yields information on potential bulk motion 
of ^^Al in the Galaxy. The ^^Al line centroid energy from 
the whole inner Galaxy is determined at 1808.92 ± 0.06 
keV from the Gaussian fit and 1809.03 ± 0.08 keV from the 
instrumental-response-convolved fit, even after considering 
uncertainties, the measured centroid energy is higher than 



the laboratory value. The effect of this blueshift for the 
whole inner Galaxy will be discussed below. 



In our determination of ^®A1 line properties, we repre- 
sent the spectra derived from detailed model fitting by a line 
component for ^^Al plus a linear component which could 
capture any underlying systematics in our background or 
sky modellings. We expect that our adjacent-energy back- 
ground model component will eliminate Galactic contin- 
uum emission, at least to first order, if the spectral shape is 
rather flat across our 40 keV energy range (1785-1826 keV). 
The diffuse gamma-ray continuum in the inner Galaxy is 
~ 2 X 10~^ ph cm~^ s~^rad~^ keV~^ in the energy band of 
1-2 MeV (Strong et al. 1999). Indeed, we find that offsets 
above zero in our spectra are rather small and negligible, 
supporting this property of our background model. 
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Fig. 8. ^^Al spectra of four segments along the Galactic plane ( —60° < I < 60°). 



4. ^^Al emission along Galactic longitudes 

The 9-year COMPTEL imaging of ■^^Al hne emission had 
already suggested some asymmetry in the inner Galaxy: 
the fourth Galactic quadrant appears somewhat brighter 
than the first quadrant (Pliischke (2000) finds a significance 
of 2.5(7 for a brightness difference). COMPTEL could pro- 
vide the image details of ^^Al in the Galaxy, but no signif- 
icant spectral information due to its spectral resolution of 
about 150 keV near the ^^Al line. SPI with its Ge detectors 
features sufRciently-high spectral resolution to allow astro- 
physical constraints from ^^Al line shapes, averaged over 
the Galaxy as discussed above, but also for different regions 
along the Galactic plane due to its imaging properties as a 
coded-mask telescope. 

In this section, we will proceed towards increasing spa- 
tial resolution along the Galactic plane, starting out from 
testing Galactic asymmetries between the first and fourth 
quadrant. ^^Al line parameters toward the different direc- 
tions of the Galactic plane are determined using separate 
sky maps covering each sky region, simultaneously fitting 
these together with our background model to the entire 
sky survey database. ^^Al line fluxes, centroid energies, and 
line widths then are derived by a simple Gaussian fit to 
the ^^Al line in the resulting spectra, as we are interested 
in relative changes between different portions of the sky. 
This will allow us to identify line shifts from bulk mo- 
tion such as expected from large-scale Galactic rotation, 
and hints for additional line broadenings in particular re- 
gions, which would refiect increased ^^Al velocities in such 
regions. A homogenous disk model (see Figure 2 right, 
-60° < I < 60°, -10° < b < 10°, scale height 200 pc 
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Fig. 9. ^^Al spectrum from the Galactic center (—5° < I < 
5°, —10° < 5 < 10°). The determined line centroid energy 
(1808.58±0.27 keV) is consistent with the laboratory value. 



) is used here to avoid a bias of sky distribution models 
along the Galactic plane, for such relative comparison. 

26A1 spectra for the 1st (0° < ^ < 60°) and 4th quadrant 
(—60° < / < 0° ) are presented in Figure 7. In the 4th quad- 
rant we note a blueshift of 0.49 ± 0.07 keV relative to the 
centroid energy of ^^Al line in the laboratory, but no signif- 
icant redshift in the 1st quadrant is apparent 0.04 ±0.10 
keV). Both spectra have width values compatible with no 
significant ^^Al line broadenings. The indicated ^^Al asym- 
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metry between the two inner Galactic quadrants appears 
again, with a flux ratio of ^ 1.3 ± 0.2. 

We proceed further towards more confined Galactic re- 
gions using four sub-maps with 30 degree width along 
Galactic longitude, to obtain ^^Al spectra for these regions 
shown in Figure 8: (1) 0° < / < 30°, (2) 30° < I < 60°, 
(3) -30° < / < 0°, and (4) -60° < I < -30°. Centroid 
energy shifts of the ^^Al line are found, -f0.15 ± 0.12 keV 
and —0.61 ± 0.09 keV in regions (1) and (3) respectively, 
as expected from large-scale Galactic rotation. The inner 
region (—30° < / < 30°) is seen to be much brighter than 
the two outer regions of the Galaxy, consistent with the 
COMPTEL -^^Al map (note that here we use a homoge- 
neous sky distribution model as the model-fitting prior). 
The indicated ^^Al emission asymmetry for the 1st and 4th 
quadrants also shows up again between regions (3) and (1), 
with a flux ratio of ~ 1.15 ± 0.18. Even further out, re- 
gion (4) may also be brighter than (2) by ~ 1.6 ±0.6. Since 
smaller regions include less ^^Al signal, these differences are, 
however, not significant. 

Challenging the imaging capability of SPI for diffuse 
and extended emission, we refine spatial structure even 
more towards smaller longitude intervals. In Figure 9, we 
show a spectrum representing the Galactic center region 
(-5° < ; < 5°, -10° < 5 < 10°), where the ^^Al line is stiU 
significant (> 4tT). The line centroid energy is determined 
at 1808.58 ±0.27 keV, and is consistent with the laboratory 
value - no line shift from bulk motion is indicated. 

Then we also derive spectra for six smaller longitude in- 
tervals of 20 degree width along the Galactic plane (—60° < 
Z < 60°, see Figure 10). The ^^Al line is stiU detected for 
the inner Galaxy (-40° < I < 40°, > 6a for each 20° 
bin region), but only marginal for the two outer regions 
(40° < / < 60° and -60° < I < -40°, < Aa). This may be 
attributed to both less ^^Al brightness and to less exposure 
in these regions, compared to the inner Galaxy. 

For the four inner regions, ^^Al line centroid shifts are 
observed. When we use the ^^Al line centroid energy of 
1808.65 keV from the laboratory as a reference (rather than 
our own determination. Fig. 9), we obtain ^^Al line centroid 
energy shifts along Galactic longitudes as shown in Figure 
11. Evidently, towards positive longitudes of ^ 0° — 40°, 
the redshift in energy is minor, ~ 0.1 keV, while for nega- 
tive longitudes, the blueshift is substantial, 0.4 keV for 
-20° < / < 0° and up to 0.8 keV for -40° < I < -20°. 
This asymmetry of ^^Al line energy shifts along the Galactic 
plane is inconsistent with simple azimuthally-symmetric 
Galactic rotation; it is clear evidence for pronounced spiral- 
arm structure in the inner Galaxy, and possibly pecu- 
liar bulk motion along the Galaxy's bar (e.g. Englmaier 
& Gerhard 1999, Hammersley et al. 2007 and references 
therein). 

Figure 11 shows the intensity distribution of ■^^Al emis- 
sion (from the same 20° regions. Figure 10) along the 
Galactic plane, adding the (longitude-range normalized) 
Galactic-Center region ^^Al intensity < 5°, Figure 9) 
for comparison. The variability of ^^Al intensity along the 
Galactic plane again is evident. 

We note that the ^"^Al fine in the region of 20° < / < 40° 
appears somewhat broadened, with a Gaussian width of 
FWHM~ 4.15 ± 0.75 keV (also see Figure 11). This may 
hint towards a peculiar ^^Al source region towards this di- 
rection, which could be associated with the Aquila region 
(Rice et al. 2006). Broadening could result from higher tur- 
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Fig. 11. Top: ^^Al line energy shifts along the Galactic 
plane ( -60° < I < 60°, with longitude bin widths of 20°, 
^^Al spectra from Figure 10) relative to the line centroid of 
the 2^A1 line (fitted energy 1808.65 keV). Center: ^^Al in- 
tensity distribution along the Galactic plane (from Figures 
9, 10). For comparison, the COMPTEL-derived ^^Al in- 
tensity profile is shown (solid line, and dashed lines when 
integrated over the same longitude bins ). Bottom: '^^Al 
FWHM (Gaussian fitting) variation along the Galactic lon- 
gitudes, a broad ^®A1 line feature is detected toward the 
longitudes 20° < ^ < 40°. 

bulence if the -^^Al source region is younger than average 
and dominated by the ^®A1 ejection from more massive stars 
(see Knodlseder et al. 2004). Further studies would be in- 
teresting, and have the potential to identify star formation 
otherwise occulted by foreground molecular clouds. 

5. Latitudinal variations of ^^Al emission 

The interpretation of ^®A1 imaging and spectral results re- 
lies on (uncertain) distances of ^®A1 sources. Along the line- 
of-sight, the detected ^^Al signal could originate from local 
star- formation complexes ('^ 100 pc), or from the nearest 
part of the Sagittarius-Carina arm (1 — 2 kpc), or from the 
Galactic center region 8 kpc), or even from the distant 
side of the Galaxy (> 10 kpc). 
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Fig. 10. ^^Al spectra of six segments along the Galactic plane ( —60° < I < 60°). Small longitude degree bin (20°) makes 
the detections of ^^Al not significant in the regions of 40° < |^| < 60°. 



In this section, we explore a possibility to resolve the 
^^Al signals for local complexes from the large scales of 
the Galaxy through different latitudinal signatures. For 
such analysis, we split a homogenous-disk model (Figure 
3 right) into sub-maps along Galactic latitudes, and de- 
rive separate ^^Al spectra for the different latitude ranges 
simultaneously through model fitting to our observations. 
We study the ^^Al emission for three intervals along lati- 
tudes, -5° < 6 < 5° (low latitudes), 5° < 6 < 20° and 
-20° < 5 < -5° (intermediate latitudes). ^'^Al emission for 
low latitudes (|6| < 5°) should be dominated by a large- 
scale origin in the Galactic disk, while ^®A1 emission at 
higher latitudes (|5| > 5°) could originate from more lo- 
cal star-formation systems such as OB associations in the 
Gould Belt. This definition is similar to the one used in pul- 
sar population studies (Wang et al. 2005). The Gould Belt 
is an ellipsoidal shaped ring-like structure delineated by the 
groups of nearby stellar groups within ~ 1 kpc, with semi- 
major and minor axes of 500 pc and 340 pc, respectively 
(Perrot & Grenier 2003). The center of this structure is lo- 
cated towards I — 130° , displaced from the Sun's location 
by about 200 pc (Guillout et al. 1998). The Vela region is 
located near the outer boundary of the Gould Belt towards 
I ~ —90° . The nearby Sco-Cen region at about 140 pc dis- 
tance towards I ~ —10°, b ~ 10° probably also belongs to 
the Gould Belt. The origin of this structure is debated, be- 
tween triggered sequential star formation propagating out- 
wards from its center, and an external triggering event such 
as a high- velocity cloud falling through the Galaxy's plane 
(see Perrot & Grenier 2003 and Poppcl 1997). 

Figure 12 displays the ^^Al intensity distribution along 
Galactic latitudes for |Z| < 60°. No ^®A1 signal is detected at 
negative latitudes (—20° < b < —5°), while weak ^^Al emis- 
sion is detected at positive latitudes (~ 2a, 5° < b < 20°). 
We also derive latitudinally-separated ^^Al emission for the 
1st and 4th quadrants separately again through model fit- 
ting, also shown in Figure 11. No ^^Al signal is detected 
for off-plane regions of the 1st quadrant, while in the 4th 



quadrant, ^^Al emission is still clearly detected (3(7) in the 
latitude region 5° < 6 < 20°. This ^^Al emission towards 
6 > 5°, I < 0° could be attributed to the nearby Sco-Cen 
star-formation complex at ~ 140 pc distance. 

We also determine the scale height of the Galactic plane 
in ^®A1 emission by comparing the fit quality for sets of two 
different plausible geometrical models for the ^^Al source 
density distribution in the Galaxy, varying their scale 
height parameter (Fig. 12 right). We use a Galactocentric 
double-exponential disk model (ExpDisk) described by 

piR,z) oc e Uo + =o (2) 

where R is the Galactocentric distance within the plane, 
and z is the height above the plane, with Rq — 3.5 kpc and 
zq as the scale radius and height parameter. Alternatively, 
we use a model which includes the spiral-arm structure 
of the Galaxy as derived by Taylor & Cordes (1993); we 
use only their components for the inner Galaxy and the 
spiral arms, with identical scale height parameter (in this 
model, the density perpendicular to the disk is described 
as sech(z/zo)). In order to avoid a bias from bright spe- 
cial regions such as Cygnus/Vela/Carina, we restrict this 
analysis to data within |/| < 60° and |&| < 30°. Fig. 12 
(right) shows the variation of log-likelihood values with dif- 
ferent scale heights for both model types. Here, the values 
for the exponential-disk models have been shifted by -1-16, 
as the spiral-arm model systematically provides a better 
description of our data. With — 21ogL being asymptotically 
distributed, we derive a scale height of 130j;yQ° pc {la) 
for the ^^Al emission in the inner Galactic disk, from the 
spiral arm model constraints. This confirms previous such 
studies based on COMPTEL data (Diehl et al. 1998). 

6. Summary and Conclusions 

The study of the ^^Al line and its details in different regions 
of the Galaxy is one of the main goals of the SPI spectrome- 
ter on INTEGRAL. Using five years of SPI data, we find the 
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Fig. 12. Left: ^®A1 intensities versus Galactic latitudes (solid lines, |/| < 60°). ^''Al emission from low latitudes (the 
Galactic disk) dominates < 5°). For comparison, ^®A1 intensities for each Galactic quadrant are also shown. Significant 
^^Al emission is detected towards 6 > 5°, I < 0°. Right: Determination of the Galactic-disk scale height. Shown is the 
logarithmic likelihood-ratio for the exponential-disk (dashed line) and spiral-arm structure models (solid line) for different 
scale height values. Log-likelihood values for the exponential-disk model have been offset by -|-16. 



26A1 signal from the inner Galaxy {\l\ < 30°, |6| < 10°) with 
a high significance of ^ 28<t. The ^^Al fiux integrated over 
the inner Galaxy of (2.9 ± 0.2) x 10"'' ph cm'^ s'^ rad'^ 
is consistent with our own and other earlier measurements, 
though significantly lower than the ones from measure- 
ments with instruments which have even more modest spa- 
tial resolution on the sky (see Figure 2 of Diehl et al. 2004). 
Taking the distance of the Sun to the Galactic center as 
i?o = 8.5 kpc, we convert the measured ^^Al fiux to a 
Galactic ^^Al mass of (2.7±0.7) Mq . We compared different 
plausible models for the sky distribution of ^^Al emission 
and how it may affect the global ^^Al flux measurement. We 
find that different sky distribution models do not substan- 
tially affect ^^Al intensity and line shapes for the integrated 
inner Galaxy, and we have used the observed variations 
of ^^Al mass values to estimate a systematic uncertainty, 
added to the statistical uncertainty from the model fitting, 
to yield the ±0.7 Mq uncertainty quoted. 

The ^®A1 line centroid energy appears blue-shifted rel- 
ative to the laboratory value of 1808.65 ± 0.07 keV if in- 
tegrated over the inner Galaxy (1809.0 ± 0.1 keV). With 
refined spatial resolution this turns out to be mostly due 
to asymmetric bulk motion which we find along the plane 
of the Galaxy, and attribute to asymmetries in inner spiral 
arms and to the Galaxy's bar. In particular, for the cen- 
tral longitude bin towards the Galactic Center, the signal 
is strong enough to allow for a rather small longitude range 
integration, and we find a ^"^"Al line centroid consistent with 
the laboratory value and thus with the absence of bulk mo- 
tion relative to the Sun. Also, here the ^^Al line appears 
as most-narrow with an upper limit of 1.3 keV {2a). The 
measured line width of ^^Al from the large-scale integrated 
inner Galaxy with an upper limit of 1.3 keV (2(t) is con- 
sistent with expectations from both Galactic rotation and 
modest interstellar-medium turbulence around the sources 
of ^^Al (turbulent velocities constrained below 160 km s~^ 
even when disregarding the effects of galactic rotation) . Our 
line width results are consistent with previous reports by 
HEAO-C (Mahoney et al. 1984), RHESSI (Smith 2003), 
but the very broad line with a width ~ 5.4 keV reported 
by GRIS (Naya et al. 1996) is clearly ruled out by our SPI 
measurements ( 4c7 significance level ). 



From our study of ^^Al emission in spatially-restricted 
regions along the plane of the Galaxy, we find: 

— ^^Al brightness appears asymmetric for the two inner 
quadrants, with a fiux ratio of the 4th quadrant to the 
1st of ~ 1.3 ±0.2 (Figure 7). 

— The ^^Al line energy varies clearly along the Galactic 
plane (Figure 11): a minor redshift (~ 0.1 keV) for pos- 
itive longitudes, but significant blueshift (~ 0.4 — 0.8 
keV) for negative longitudes. 

— The ^^Al line towards the direction of 20° < Z < 40° 
shows a hint for additional line broadening. 

— The scale height of Galactic-plane ^^Al emission is 

pc (Icr), determined towards the inner Galaxy 
(|Z|<60°). 

— There is a strong hint for ^®A1 emission in the fourth 
Galactic quadrant at intermediate latitudes Z < 0°, b > 
5° (Fig.ll). 

This leads us to several astrophysical implications and 
the "bar" structures: 

— The disk of our Galaxy apparently is not azimuthally- 
symmetric. The ^^Al brightnesses for the two inner 
quadrants appear different, both from our INTEGRAL 
and from earlier COMPTEL results, favoring the fourth 
quadrant. ^^Al line shifts towards the blue in such 
brighter regions are more pronounced than the redshifts 
on the fainter side of the disk, unexpected from a sim- 
ple and symmetric Galaxy model and its large-scale ro- 
tation propc;rtics. These phenomena may reflect varia- 
tions in the space distribution of young stars, possibly 
related to the inner parts of the spiral arms and their 
interfaces to the "molecular ring" and the "bar" struc- 
tures. Large non-circular motions have been seen in HI 
and CO observations (Mulder & Liem 1986; Gerhard & 
Vietri 1986), and also from the NIR light distribution 
(Binney et al. 1997), in addition to source count asym- 
metries (Nikolacv & Weinberg 1997) and non-symmetric 
gas dynamics (Englmaier & Gerhard 1999). The bar 
is most clearly traced in infrared emission from dust 
(Marshall et al. 2008), which suggests a position an- 
gle of ~ 25° of the bar with the near end pointing 
towards us in the first quadrant, and a total length 
of about 4 kpc. The transition regions between spiral 
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arms and bar are likely to incur star formation (Verley 
et al. 2007), and affect the dynamics of gas and stars. 
Simply superimposing a homogeneous star-forming bar 
to the rotating spiral would, however, not agree with 
our results - the inner Galaxy structure may be more 
complex. Localized star-forming regions could lead to 
peculiar motion of hot gas ejected from winds and super- 
novae, which might dominate over large-scale Galactic 
rotation. Such activity could be the cause of the ob- 
served asymmetry of the ^^Al line intensities and line 
energy shifts in the inner Galaxy. On the other hand, 
other more nearby ^^Al source regions could be respon- 
sible for these irregularities, such as the nearest part 
of the Sagittarius-Carina arm in the fourth quadrant, 
or regions/complexes attributed to the Gould Belt such 
as the Scorpius-Centaurus-Lupus groups. Refined ^^Al 
studies and their combination with astrometry from 
other tracers of Galactic structure could help to under- 
stand the role of our Galaxy's bar. 

— The ^^Al emission asymmetry (1.3 ± 0.2) between 
fourth and first quadrant of the Galaxy is lower than 
the intensity contrast of l.Slg'g reported for positron 
annihilation emission from the disk of the Galaxy 
(Weidenspointner et al. 2008). ^^Al by itself releases 
a positron in 82% of its decays; it is uncertain, how- 
ever, how this translates into annihilation photons, from 
the variety of slowing down and annihilation processes 
which determine the fate of positrons in interstellar 
space. Both spatial and temporal variations and non- 
linearities scaling with gas density may occur. 

— The scale height of Galactic-plane ^^Al emission is sig- 
nificantly larger than the molecular-gas disk scale height 
of ~ 50 pc, yet significantly smaller than the "thick 
disk" part of the Galaxy. It is consistent with ^^Al being 
ejected from star-forming regions, and partly extending 
more towards the Galactic halo where gas pressure is 
lower than within the plane of the Galaxy ( "champagne 
flows"). 

— Localized regions may deviate in interesting detail from 
the large-scale averaged properties of ^^Al source re- 
gions. The direction of 20° < I < 40° corresponds to the 
Aquila region, and our hint for additional line broaden- 
ing may be due to increased interstellar turbulence from 
stellar-wind and supernova activity at the characteris- 
tic age of stellar groups of that region, which may have 
created a supershell of substantial size (320 x 550 pc, 
see Maciejewski et al. 1996). Similar arguments are be- 
ing explored for the Cygnus region, based on earlier 
hints of peculiar ^^Al emission (Martin et al. in prepa- 
ration). The hint for ^^Al emission at / < 0°, 6 > 5° 
(Fig. 12) may be attributed to the relatively nearby star 
forming complexes of the Sco-Cen association (de Geus 
1992). Sco-Cen and several nearby stellar groups are 
attributed to the larger structure of the "Gould Belt" , 
which is suggested to have been more actively forming 
stars during the last 30 Myr (Grenier 2000, Perrot & 
Grenier 2003). These examples indicate that spatially- 
resolved ^^Al emission properties may enable new diag- 
nostics of the interactions of massive stars with their 
surroundings, when combined with other astronomical 
constraints on cold gas and stars of such regions. 

INTEGRAL/SPI will continue to accumulate more data 
for several years, to cover more sky regions in the outer 



Galaxy and at intermediate latitudes. This will allow us 
to deepen the studies reported in this paper, and to refine 
spatial information on the ^^Al line and/or increase the sig- 
nificance of the results reported above, as signal statistics 
increases and instrumental-background models are tight- 
ened. These studies will help us to better understand the 
properties of ISM near groups of massive stars, and the 
bulk motion of gas in the inner Galaxy from Galactic ro- 
tation and other peculiar kinematics. Nearby ^^Al sources 
may be also discriminated with better exposure towards 
candidate regions, improving our determination of the ^^Al 
mass in the Galaxy and towards regions where the stel- 
lar census is known to a better degree. Studies of ^^Al from 
nearby star-formation regions (e.g., Cygnus, Vela, Sco-Cen, 
Orion) are a promising diagnostic for the massive star ori- 
gin of Galactic ■^^Al and kinematics of ^^Al ejecta in ISM. 
From this, important constraints on massive-star and su- 
pernova nucleosynthesis are obtained, refining our models 
of their complex interiors (see e.g. Woosley and Heger 2007, 
or Chieffi and Limongi 2006). 
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